Finite Temperature Momentum Distribution of a Trapped Fermi Gas 
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We present measurements of the temperature-dependent momentum distribution of a trapped Fermi gas con- 
sisting of ''"K in the BCS-BEC crossover regime. Accompanying theoretical results based upon a simple mean- 
field ground state are compared to the experimental data. Non-monotonic effects associated with temperature, 
T, arise from the competition between thermal broadening and a narrowing of the distribution induced by the 
decrease in the excitation gap A(r) with increasing T. 



PACS numbers: 03.75.Hh, ()3.75.S.s, 74.20.-z 

The recent discovery of superfluidity in trapped fermionic 
gases has paved the way for arriving at a much deeper un- 
derstanding of the phenomena of superfluidity and supercon- 
ductivity ClISISHIllSlllil. This may ultimately have 
application in high temperature superconductors S Uoll . nu- 
clear, astro- and particle physics. In these trapped gases one 
has the ability to tune the strength of the attractive interac- 
tion that leads to pairs of fermions (Cooper pairs), which 
then Bose condense. As the pairing strength is increased a 
smooth crossover from BCS behavior to Bose-Einstein con- 
densation (BEC) occurs. This tunability is a remarkable fea- 
ture of trapped atomic gases and arises from a phenomenon 
known as a Feshbach resonance. In general the resulting su- 
perfluid state is more complex than that of simple BCS the- 
ory. Similarly the normal state is expected to be quite different 
from its BCS analogue, since in general pairing takes place at 
higher temperatures (T*) than the condensation temperature 
(Tc). 

The ultracold gases have, thus, presented us with an oppor- 
tunity to investigate in a more complete fashion all aspects 
of fermionic superfluidity. Because of their charge neutrality 
and because of trap confinement effects there exists a differ- 
ent set of tools for their experimental investigation. Among 
these tools are measurements of the real space and momentum 
space distributions of the atomic fermions Illlfl2lll3lfl4lllal . 
At finite temperatures these distributions must change signifi- 
cantly in ways that reflect both pairing and possibly the onset 
of phase coherence. 

In this paper we study both experimentally and theoreti- 
cally the temperature dependence of the momentum distribu- 
tion of fermionic atoms within a trapped Fermi gas. Recent 
measurements of this momentum distribution near T — 
showed a dramatic broadening as the pairing strength in- 
creased 1 13|], as is qualitatively consistent with theoretical 
calculations I13l V\M . In the present paper we extend these 
measurements by varying the temperature from below to well 
above the theoretically predicted values of Tc\ this allows us 
to probe the momentum distribution in the normal phase as 
well as the superfluid phase. 

In general, in a T = fermionic system the physical phe- 
nomenon that controls the momentum distribution is pairing. 
In the BCS theory the momentum distribution of a homoge- 
neous sample shows a slight smearing of the Fermi surface 
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due to pairing. This effect is very small and associated with 
the gap parameter A. In the BCS-BEC crossover, however, 
one greatly increases the interaction strength that leads to pair- 
ing. It is conventional to parameterize the state within this 
crossover in terms of the dimensionless product A;^a, where a 
is the two-body s-wave scattering length between atoms and 
k''p is the Fermi wavevector in the noninteracting limit. While 
in the BCS limit (l/kpU — > — cxo) it is difficult to observe the 
slight change in the momentum distribution, as the interaction 
strength increases this broadening grows. Indeed, at unitarity 
(l/kptt = 0) the broadening is no longer small, but instead 
comparable to the Fermi energy. Finally, in the BEC limit 



(l/k°p, 



cx)) the atomic momentum distribution becomes 



extremely broad and corresponds to the square of the molecu- 
lar wavefunction in momentum space. 

At finite T one expects a similar broadening of the momen- 
tum distribution as the interaction is increased. However, the 
extent of this broadening is determined by the temperature de- 
pendence of the excitation gap A(r); it is maximum at T = 
and should disappear as A(r) goes to zero around T*, where 
all the fermion pairs are broken. This effect of pamng on the 
distribution will occur in addition to the usual thermal broad- 
ening of the momentum profile that occurs in a Fermi gas. 
While the thermal broadening increases with T, the pairing 
induced broadening decreases with T. As we will show be- 
low, these two competing effects lead to a non-trivial T de- 
pendence in the Fermi gas momentum distribution. 

While the discussion thus far has centered on the expecta- 
tion for a homogeneous Fermi gas, the experiments presented 
here take place in a harmonic trapping potential. In such a 
potential the local Fermi wavevector varies with position r, 
so that the integrated momentum distribution for a noninter- 
acting Fermi gas in a trap is significantly different from that 
for a homogeneous case 1 13]. Thus, in addition to the local 
broadening due to pairing, the trapping potential compresses 
the density profile and thereby enlarges the overall momentum 
distribution. 

To perform our experiments, we follow the techniques de- 
scribed in Refs. ITtI [lai to create an ultracold ^°K gas in 
the BCS-BEC crossover. Briefly, starting with a cold gas we 
slowly ramp a magnetic field to approach a wide Feshbach 
resonance at 202.1 ± 0.1 G |31- The rate of the magnetic 
field ramp, (6.5 ms/G)~^, is slow enough to ensure adiabatic- 



ity with respect to many-body time scales in the system. We 
then probe the momentum distribution of atoms in this final 
state using the experimental sequence developed in Ref. I13ll . 
In this sequence the atom gas must expand freely without 
any interatomic interactions; this we achieve through a fast 
magnetic-field ramp (here « 8 /is/G) to a ~ lllll . This 
is followed by the standard technique of time-of-flight expan- 
sion and optical absorption imaging. 

To obtain the data in this paper we repeat the process above 
for gases at a variety of interaction strengths and at a vari- 
ety of temperatures. To vary the temperature of our initial 
Fermi gas we increase the depth of our optical dipole trap af- 
ter evaporative cooling and heat the gas by modulating the 
optical trap strength. To characterize the temperature of the 
gas we measure the temperature compared to the Fermi tem- 
perature in the non-interacting limit, {T/Tp)'^. The physical 
temperature at a given kpU can then be extracted using the the- 
oretical thermometry described in Ref. fly]; this thermometry 
is based upon entropy conservation in the adiabatic magnetic 
field ramp. 

In Fig. 1, we plot experimental azimuthally averaged 
momentum distributions as a function of k/kp. From top 
to bottom, each panel corresponds to a fixed temperature 



(T/Tp 



0.11, 0.2, 0.3, 0.5, 0.7, and 1.0, with variable 



1/k^a. The scattering strength 1/fc^o = -71, -0.5, 0.0, and 
0.5 represent the noninteracting Fermi gas, near-BCS, unitary, 
and near-BEC cases, respectively. For each set of measured 
distributions in Fig. 1 at a particular {T/Tp}'^, parameters 
such as trap strength and initial gas density are fixed. How- 
ever, these parameters vary among the different sets at con- 
stant (T/Tp)'^ (panels in Fig.Q. For example, for the data 
at [T/Tp)'^ = 0.11 the peak density, for atoms in one spin 
state, in the weakly interacting regime is n^^ = 1.4 x 10^"^ 
^^ and E^ = 0.56 /iK. For the hottest data n^j. decreases 



cm 



to 6 X 10'^ cm-'' and E'^, = 0.79 /^K. 

On top of the data points in Fig. 1 we show a fit to an empir- 
ical function, which we utilize both for normalization and for 
later data analysis. We apply the following two-dimensional 
(2D) surface fit to the measured optical depth (OD) of the ex- 
panded Fermi gas. 



" /.92(-C), (1) 



ODix,y) = ODpk92[-Ce 



=>) Jo -^T3i?7J = ^'=1 ^ ^' *^ P°^y" 
logarithmic function, ( is the fugacity, a-^ are proportional to 
the Fermi gas temperature and related to the expansion time t. 
This form can be derived microscopically in the limit of weak 
interactions, where ( — > e^/^^^ . While this equation is only 
physically valid for an ideal Fermi gas, empirically we find 
that it fits reasonably well to data throughout the crossover 
llol . In this way the 2D momentum distribution is given 
by N2D{kx,ky) = A OD{hkxt/m,hkyt/m). Here A is a 
normalization constant such that J ."^J:.^ N2D{kx, ky) = Na, 
where Na is the total number of atoms. 

The experimental results of Fig. 1 show that, at {T/Tp}'^ = 
0.11, as the system passes from an ideal Fermi gas to near- 
BEC the momentum distribution widens significantly, just as 




FIG. 1: (color online) Evolution of the momentum distribution 
N2D{k) with the interaction strength \/k%a from noninteracting to 
near-BEC. The two columns compare experiment (left) and theory 
(right). Different rows correspond to different values of {T/Tf)^. 
The lines in the experimental plots correspond to a fit to Eq.Q The 
value of k% used to create the experimental plots was determined 
through experimental measurements of the particle number and trap 
strength, which lead to systematic uncertainties in fc^ of up to 10% 
for these data. While the raw data are shown here, the normalization 
applied in Fig.|3|removes this systematic error 



the T = case shown in Ref. |131. As {T/Tp)° increases the 
effect of pairing becomes less pronounced, and at sufficiently 
high T the curves coalesce. Here the gap A(r) is small com- 
pared to T, and the system basically behaves as a classical gas 
of atoms. 

Theoretical calculations for the momentum distributions in 
the crossover are presented in the right column of Fig. [2 
The calculations are based on a generalized mean-field the- 
ory 1 9, 20], which is consistent with the BCS-Leggett ground 
state |21|. While T — Q momentum profiles of this state have 
been discussed in the literature 1 13, 22], here we include fi- 
nite temperature effects. An important aspect of this theory is 
that the fermionic excitation gap A(T) becomes distinct from 



the superfluid order parameter A^c (T) at finite T (except in 
the strict BCS regime). This point is illustrated in Fig. EJa). 
Here the trap-averaged value of the excitation gap, {A^y^^, is 
shown as a function of {T/Tp}'^ for the near-BEC [l/kpa = 
0.5), unitary (l/kpa = 0), and near-BCS {l/kpa = -0.5) 
cases from top to bottom. The arrows indicate the calculated 
value of T°. We find that A(T) has a gentle onset at high T, 
corresponding to the pair formation temperature T*, which is 
generally much higher than Tc- 

The momentum and density distribution, along with the 
chemical potential /i, the gap parameter, and the order param- 
eter, are self-consistently determined using the local density 
approximation. The local momentum distribution is given by 



nk{r)^2[vl{l-f{E^))+iLlf{E^)], 



(2) 



where ul,vl = [1 ± (ek - fJ.{r))/Ek{r)]/2 are the BCS 
coherence factors; ek ~ h^k'^/2m is the kinetic energy of 
free fermions, _Ek — ^(ek — l^ii~)y + ^^ the fermionic 
quasiparticle dispersion, and f{x) the Fermi-Dirac distribu- 
tion function. Then the trap-integrated momentum distribu- 
tion N{k) and its 2D projection are given, respectively, by 



iV(fc)= d^rnkir) and N2D{k) 



dkz 

2^ 



N{k). (3) 



Note, it is the 2D distribution A^2u(fc) that is directly compa- 
rable with our experimentally measured optical depth. It can 
be seen that the effects of superfluidity and of pairing more 
generally enter only through Ek which, in turn, depends upon 
A(T). Thus the momentum distribution represents, in effect, 
the behavior of the fermionic excitation gap, not the superfluid 
order parameter Age- Details of the theoretical formalism can 
be found in Ref . ^9J . 

Near T = it has been shown that the simple mean-field 
ground state we use here semi-quantitatively describes the 
momentum distribution and effective kinetic energy I 13ll23ll . 
Thus, this ground state is an excellent starting point for fi- 
nite T comparisons between theory and experiment. Quanti- 
tatively, weaknesses of the mean-field ground state in repro- 
ducing the momentum distribution exist, as seen in Ref. I23ll . 
These weaknesses include that, in the BEC regime, the inter- 
boson scattering length is overestimated by roughly a factor of 
three. However, this does not substantially affect the momen- 
tum distribution of atoms we address, in part because for the 
kpa considered here we confine our attention to the fermionic 
regime (/i > 0). (For l/k^a — 0.5, our calculation shows 
fj,{r = Q)/E% = 0.23 at Tc.) Secondly, the theory does not 
include the Hartree self energy. We also do not take into ac- 
count the finite rate of the magnetic-field ramp above the Fes- 
hbach resonance as was done theoretically in Ref. 1 13|. This 
finite ramp rate leads to slight redistribution of the momentum 
so that high energy spectral weight is somewhat suppressed. 
This effect is crucial for accurately comparing to the effec- 
tive kinetic energy 11131 l23ll . but for our ramp rates and in the 
fermionic regime we expect the effect on our distribution com- 
parisons to be small. 

We now focus on the comparison between our mean-field 
calculations and the experiment shown in Fig. 1. Overall we 
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FIG. 2: (color online) (a) Temperature dependence of the average 
excitation gap (A^)^" at 1/kpa = 0.5 (top blue curve), (red) 
and -0.5 (bottom black curve). The solid circles indicate where this 
effective gap coincides with the temperature. The arrows indicate the 
transition temperatures T^. (b) Theoretical momentum distributions 
N2D{k) at a variety of temperatures for l/kpa = 0.5. 



find semi-quantitative agreement. However, there are discrep- 
ancies that can be seen in the profiles especially in the middle 
range of 1/kpa. This can be attributed mainly to the neglect 
of the Hartree term in the theoretical formalism that was noted 
above. If we were to include the Hartree correlation our calcu- 
lation of the chemical potential /i would decrease; this would 
result in a wider spread between the noninteracting and inter- 
acting profiles in the theoretical portion of Fig.^and, thus, a 
better quantitative agreement with experiment. 

We now turn to understanding the T dependence of the dis- 
tribution at a fixed value of l/k^a. In this analysis one ex- 
pects a competition between conventional thermal broadening 
of the momentum profiles and narrowing of the profiles due 
to the decrease in size of A(r) with increasing temperature. 
To illustrate this quantitatively we turn to Fig. OJa). The solid 
circles in this figure indicate where the effective gap coincides 
with the temperature. This corresponds roughly to the tem- 
perature at which the broadening as a function of T should 
display an inflection point. 

In Fig. |3b) we plot the theoretical results for the momen- 
tum distribution at a fixed l/kpU ~ 0.5. Indeed, there is 
a change in the temperature dependence of the curves which 
occurs roughly when (A^)^'^ = ksT. Below this tempera- 
ture the behavior of the distribution tends to be dominated by 
that of the excitation gap, and the distribution initially narrows 
with increasing T. Above this temperature, thermal effects 
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FIG. 3: (color online) Comparison of the half width at half height 
(HW) as a function of (T/Tf)''' for different values of 1/kpa from 
noninteracting (bottom curve) to near-BEC (top curve) between (a) 
experiment and (b) theory. The HW in the experimental case is 
normalized to eliminate dependence upon the experimental determi- 
nation of Ep (see text); this results in perfect agreement with theory 
in the non-interacting case (black line). 



dominate and the distribution widens with increasing T. We 
find analogous non-monotonic behavior at the unitary limit. 

To test this theoretical prediction of a non-monotonic tem- 
perature dependence of the broadening, we plot the half width 
at half height (HW) of the distributions as a function of 
{T/TpY' . The experimental HW is determined from a fit of 
the data to Eq. fusing the curves in Fig. n]as well as addi- 
tional similar data. For this analysis we have eliminated the 
dominant uncertainty associated with the determination of E^ 
(and hence /c^) by applying a multiplicative correction factor; 
this factor is simply the ratio of the calculated and measured 
half widths for the weakly interacting limit. This is done for 



each set of data taken at fixed {T/Tf)°. 

Figure |3a) shows the experimental results; the different 
symbols represent groups of data with an average value of 
l/kpG indicated by the legend. The black line represents 
the theoretical dependence of the HW in the non-interacting 
(a = 0) limit. As expected for this limit, in the classical 
regime the width scales linearly with the temperature; the 
width then levels out to a finite value due to Pauli pressure. As 
the interaction is increased (black circles to red stars) the de- 
crease of the width with decreasing T becomes less dramatic, 
and at the lowest T the width actually begins to increase as a 
function of T. This temperature non-monotonicity is particu- 
larly apparent on the BEC side of the Feshbach resonance. 

The experimental result of Fig. [Sja) compares favorably 
with its theoretical counterpart (Fig. |3lb)) in that they both 
display non-monotonic T dependence. Note that the theoret- 
ical HW at low T and weak interactions is smaller than its 
experimental counterpart; in this regime the width drops be- 
low the noninteracting curve due to the absence of the Hartree 
term and changes in the shape of the distribution. This as- 
pect of the theory result again accentuates that, while we have 
found the simple mean-field theory to be an excellent starting 
point for finite temperature studies, future quantitative studies 
will require theories that incorporate the Hartree term. 

In summary, we have studied the momentum distribution 
of trapped fermionic atoms at finite temperature by a detailed 
comparison of theory and experiment. Our results show that 
there is a competition between the temperature dependence of 
the fermionic excitation gap and thermal broadening, which 
leads to non-monotonicities in the temperature dependence of 
the momentum profiles. Since temperature is often difficult 
to determine experimentally, systematic studies at non-zero T 
are just beginning in these ultracold gases. Here, by working 
with a theoretical temperature scale (set by the entropy), under 
conditions of an adiabatic sweep, we are able to show semi- 
quantitative agreement between theory and experiment using 
a simple mean-field theory. 
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